Reduced fidelity in Kitaev honeycomb model 
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We study reduced fidelity and reduced fidelity susceptibility in the Kitaev honeycomb model. It 
is shown that the nearest-two-site reduced fidelity susceptibility manifest itself as a peak at the 
quantum phase transition point, although the one-site reduced fidelity susceptibility vanishes. Our 
results directly reveal that the reduced fidelity susceptibility can be used to characterize the quantum 
phase transition in the Kitaev honeycomb model, which suggests that, despite its local nature, the 
reduced fidelity susceptibility is an accurate marker of the topological phase transition when it is 
properly chosen. 
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I. INTRODUCTION 

The quantum phase transition (QPT), which is a phase 
transition driven purely by quantum fluctuations and 
which occurs at zero temperature, is believed to be an 
important concept in condensed matter physics [l|, |2j. It 
is defined as the nonanalytic behavior of the ground- 
state properties, and thus reflects quantum fluctuations, 
which differentiates it from the temperature-driven ther- 
mal phase transition that reflects thermal fluctuations. 
Most QPTs can be described within the traditional 
symmetry-breaking formalism; however, there are also 
exceptions which can only be characterized by topolog- 
ical order [3j. In these topological quantum phase tran- 
sitions (TQPTs), the local perturbation effects may be 
exponentially suppressed, as was observed in fractional 
quantum Hall phases [3J. 

This QPT has appeared in a number of unrealistic 
models and is also related to many realistic systems such 
as high temperature superconductors |J|. Despite many 
theoretical examples that show the existence of the QPT, 
there is still no definite way to mark it in terms of a 
local-order parameter. Recently, much attention J5H18| 
has been drawn to the use of fidelity to mark the QPT ( 
for a review, see Ref. [19| ). For example, the fidelity ap- 
proach to the TPQT occurring in the Kitaev toric model 
have been used in Refs. [3 [ff, and the fid elity between 
two states in the Kitaev honeycomb model[2(| has been 
studied in Refs. @, Q. E. Eriksson and H. Johannesson 
proposed that several TQPTs were accurately signaled by 
a singularity in the second-order derivative of the reduced 
fidelity [9]. Fidelity is suggested to be particularly suited 
for revealing a QPT in the Bose-Hubbard model[10|. 
Moreover, the ground-state fidelity and various correla- 
tions to gauge the competition between different orders 
in strong correlated system were also discussed (lljj . 

Fidelity is a concept borrowed from quantum- 
information theory and is defined as a measure of simi- 
larity between two quantum states. Intuitively, it would 
be a good marker, since the structure of the ground state 
should suffer a dramatic change through the phase tran- 



sition point. As fidelity is purely a quantum-information 
concept, the advantage of using fidelity is that no a priori 
knowledge of any order parameter and changes of sym- 
metry of the system need to be assumed. This idea has 
been used to analyze the QPT, and its validity has been 
confirmed in many systems [19]. 

Meanwhile, the role of the leading term of fidelity was 
explored [13, [l4j . The fidelity susceptibility, which is the 
second derivative of the fidelity with respect to the driv- 
ing parameter, was introduced to study the QPT[T4|. It 
has been pointed out that the fidelity susceptibility is 
actually equivalent to the structure factor (fluctuation) 
of the driving term in the Hamiltonian on the perturba- 
tion level. In this case, if the temperature is chosen as 
the driving parameter of thermal phase transitions, the 
fidelity susceptibility, which is extracted from the mixed- 
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state fidelity between two thermal states 21], is simply 
the specific heat. From this point of view, using the fi- 
delity to deal with the QPT seems still to be within the 
framework of the correlation function approach, which 
is intrinsically related to the local order parameter. Ex- 
cept for the global fidelity, which measures the differ- 
ence between ground states differing slightly in the driv- 
ing parameter, the reduced fidelity, which describes the 
difference of the mix-state of only a local region of the 
system of interest, was also introduced. It has been 
demonstrated that the reduced fidelity should be useful 
in many ordinary symmetry-breaking QPTs as well as 
some TQPTs jl|2l|2i|. The reduced fidelity reveals in- 
formation about a change in the inner structure for a sys- 
tem undergoes a QPT, and it is significant to investigate 
the behavior of the reduced fidelity susceptibility in both 
critical and noncritical regions. Recently, it was shown 
that the TQPT in the Kitaev spin model can be char- 
acterized by nonlocal-string order parameters [23I425J . so 
it would be interesting to see whether it is able to mark 
this TQPT with a local quantity; namely, the reduced 
fidelity. 

In this work, we study the reduced fidelity and re- 
duced fidelity susceptibility in the Kitaev honeycomb 
model. Because of its potential application in topolog- 
ical quantum computation, such a model has been the 
focus of research in recent years, despite its being rather 
artificial[H HBUlj. We show that the reduced fidelity 
susceptibility of two nearest sites manifest itself as a peak 



at the QPT point, although the one-site reduced fidelity 
susceptibility vanishes. Our results directly reveal that 
the reduced fidelity susceptibility can be used to char- 
acterize the QPT in the Kitaev honeycomb model, and 
thus suggest that the reduced fidelity susceptibility is an 
accurate marker of the TQPT when it is properly chosen, 
despite its local nature. 

This article is organized as follows. The model and the 
calculation of reduced fidelity are presented in Sec. II. 
The results of the reduced fidelity and fidelity suscepti- 
bility are discussed in Sec. III. A summary is given in 
Sec. IV. 



II. REDUCED FIDELITY AND FIDELITY 
SUSCEPTIBILITY 



The two-site reduced density matrix can be expressed as 

m 

p( ri ,r 2 ) = Tr' (\g)(g\) 



1 3 

a.a'—O 



9)Oi*2, (7) 



where a a (a a ) are Pauli matrices a x , a v , and a z for a 
(a')= 1 to 3, and the unit matrix for a (a')= 0. 

Now we consider the nearest-two-site reduced density 
matrix, (i.e., the reduced density matrix of an x-bond. 
It can be proved that all but two parts of the density 
matrix are nonzero: 

p(r, r + x)=- (g \a*a* +x \ g) a*a* +x + -I r I I+ *. (8) 



The Kitaev honeycomb model was first introduced by 
Kitaev to study the topological order and anyonic statis- 
tics, which describe a honeycomb lattice system with a 
spin 1/2 located at each site. The Hamiltonian consists 
of an anisotropic nearest-site interaction with three types 
of bonds: J x , J yi and J z . The model Hamiltonian is 



x-bonds 



n = -J, V oiol~J, V M-J, V rfoi 



y-bonds 



2-bonds 



Jxflx Jy-tly J z ±l z , 



(i) 



where j and k denote the two ends of the correspond- 
ing bond, and J a and a a (a = x,y,z) are dimensionless 
coupling constants and Pauli matrices respectively. 

Bellow, we solve this model within the method de- 
veloped by Kitaev[20]. Here, we would like to mention 
that Z. Nussinov and G. Ortiz have also illustrated how 
Kitaev's honeycomb model may be solved directly |29|. 
Within this approach, the above Hamiltonian can be 
rewritten using fermionic operators as 



H 
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A 2 
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K* 
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(2) 



with 



If we assume the system has translational symmetry, we 
can calculate the average 

1 1 ft T? 

{ 9 \a x r a: +x \g) ^ - (g\H x \g) = --^. (9) 



V- 



L 2 dJ x 



Combining Eqs. (j5]) and ©, the reduced density matrix 
can be expressed as 

C x = (s|ofof+x|s) 
1 dE 



L 2 8J X 



ly^fii d £ g A q 

7\2 Z^ I E7 f> T h R_. 



A a 8A C 



L 2Z ^\Eq dJ x Eq dJ x 



= -y 

r,2 Z^ 



e q ■ cos q x + A q • sin q x 
E~ a 



(10) 



Similarly, for the y and z bonds, we can also obtain 
p(r, r + y) = -(g \a y r a y r+y \g) cr^cr^ +y + -I r I r+y , 
p(r,r + z) = - (g \a z r a z +z \g) a z a z ^ + -I r I 



r+z^ ,- t r-'r+z, 



(11) 



where 



£ q = Jx cos q x + J y cos q v + J Zl 
A q = J x sin q x + J y sin q y . 



The momenta take the values 



q x ( y ) 



2nn 

~T ,n ~ 2~ 



L-l L-l 



J 1 



(3) 



(4) 



where L is an odd integer; then the system size is N 
2L 2 . 

Therefore, we have the ground state 



i^^n^Lio), 



with the ground-state energy 



s ° = -EaM+ a I = -E^ 



(5) 



(6) 



C y = (g\ay<T y r+y \g) 
1 dE Q 



L 2 dJ y 
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L 2 
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E n 



C z = (g\a*a* +z \g) 
1 dE 



L 2 dJ z 



= -y 



^ 2 ?Uqj 



£q d€g Aq (9A q 

Eq 8J Z Eq dJ z 



(12) 



After arriving at a diagonal reduced density matrix, 
it is easy and straightforward to calculate the reduced 

fidelity 



1 .00000 



F = Y,^P^i 



(13) 



Thus we have the expression 



1, 



Fa = -[(l + C a )(l + C' a ) + (l-C a )(l-C> a )},(U) 

where C a are C x , C y , and C z for x, y, and z bond, re- 
spectively, and F a are F x , F y , and F z . 

Now we calculate the reduced fidelity susceptibility \f- 
For diagonal density matrices, the fidelity susceptibility 

is 



V^ (dj aPii ) 2 ^(d Js C a ) 2 



^Pi 



^Pi 



(15) 



where J s is the driving parameter. 

With this expression, we first set J x = J y and select 
J z as the driving parameter, and then we have 



I -J, 



■ cos q x 



\-J z 



cos q y + J z 



de c 



1 1 

^ q " dJ ~2 C ° Sqx ~ 2 COBQv+ ' 

1 - Jz . 1 - Jz . 
^ q = — 7, — smq x -\ — smq y , 



q dJ z 



-smq x - -sings,, 



q ■■ dJ z ~ E^ • £ q ~ q ' 



£ i. e ' + ^ A ' 



(16) 



and 

dj,C x = 72 E ^2 



q q 



[(£ q e q - £ q e q ) cosfe + (£ q A q - £ q A q ) sin fe ] 



q q 



[(-Eq< - S q e q ) cos q„ + (£ q A q - £ q A q ) sin g„] 

^■r ■ = J^^ Mz^fi , ( i7) 



L 2 



E 2 

^q 



Also, we can set J z = 1/3 and make J x the driving 
parameter, and the formula can be achieved in a similar 
way. Then, we can study the reduced fidelity and reduced 
fidelity susceptibility depending on which parameters we 
are interested in. 



III. RESULTS AND DISCUSSIONS 

To study the reduced fidelity susceptibility of the 
Kitaev honeycomb model, the parameters shall be re- 
stricted to the J x + J y + J z = 1 plane. According to 
Kitaev [20], this plane is divided into a gapped phase and 
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FIG. 1: (Color online) Two-site reduced fidelity as a function 
of J z along the dashed line shown in the triangle for the x, y 
and z bonds. 
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FIG. 2: (Color online) Two-site reduced fidelity as a function 
of J x along the dashed line shown in the triangle for the x, y, 
and 2 bonds. 



a gapless phase. These two phases are separated by the 
triangle line connecting the J x = 1/2, J y = 1/2, and 
J z = 1/2 points. As we already arrived at the expression 
of the reduced fidelity and the reduced fidelity suscepti- 
bility, now we begin to calculate these quantities numeri- 
cally. We start from the nearest-two-site reduced fidelity, 
since there are three types of nearest sites for a honey- 
comb lattice. In Figs. [1] and [2j the two-site reduced 
fidelity with system size L = 3001 for all three different 
bonds are shown along two orthogonal parameter direc- 
tions, J x = J y and J z — 1/3, which are marked in the 
inset of the figure. In Fig. [lj the parameter is chosen 
to be J x = J y , and the reduced fidelities are plotted as 
a function of J z , where we can see that the x bond and 
y bond reduced fidelities are the same, whereas the z- 
bond reduced fidelity differs. In Fig. [51 the parameter 
is chosen to be J z = 1/3 and the reduced fidelities are 
plotted as functions of J Xl where we can see that the 
x— and y— bond reduced fidelities are symmetric with 
respect to J x = 1/3, and the z— bond reduced fidelity 
differs. However, in these two cases, among all three 
bonds, the reduced fidelity manifests a dip at the phase 
transition point. It is quite clear that the two-site re- 
duced fidelity can serve as a signature for the TQPT in 
the Kitaev honeycomb model. 
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FIG. 3: (Color online) Two-site reduced fidelity susceptibility 
as a function of J z along the dashed line shown in the triangle 
for x, y, and z bond. 
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FIG. 4: (Color online) Two-site reduced fidelity susceptibility 
as a function of J x along the dashed line shown in the triangle 
for x, y, and z bond. 



In order to illustrate the ability of reduced fidelity to 
mark the TQPT in Kitaev honeycomb model, we cal- 
culate the reduced fidelity susceptibility and present the 
results of the reduced fidelity susceptibility for a system 
size L=3001 and for different bonds along J x = J y and 
J z = 1/3 in Fig. [3] and Fig. @] As expected, the reduced 
fidelity susceptibility shows a sharp peak at the TPQT 
point, which would be a clear sign of the phase transition. 
Comparing with the global fidelity susceptibility, the os- 
cillation in the B phase seems to be missing. Since it 
has been shown that the oscillating of the global fidelity 
susceptibility might be related to the long-range correla- 
tion function, this disappearance seems to represent the 
locality of the bond-reduced fidelity susceptibility. 

To further understand the properties of the reduced 
fidelity susceptibility, we study its scaling behavior. In 
Figs. 5-9, the reduced fidelity susceptibility of different 
bonds at different parameter lines with different system 
sizes L are presented. Figure 5 ( Fig. 6 ) shows the 
x— bond ( z— bond ) two-site reduced fidelity susceptibil- 
ity as a function of J z , and Fig. 8 ( Fig. 9 ) shows the 
x— bond ( z— bond ) two-site reduced fidelity susceptibil- 
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FIG. 5: (Color online) a;-bond two-site reduced fidelity sus- 
ceptibility as a function of J z along the dashed line shown in 
the triangle for the various system sizes L=301, 1001, 3001, 
6001. 
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FIG. 6: (Color online) z-bond two-site reduced fidelity sus- 
ceptibility as a function of J z along the dashed line shown in 
the triangle for the various system sizes L=301, 1001, 3001, 
6001. 



ity as a function of J x , along the dashed line shown in the 
triangle for the various system sizes L=301, 1001, 3001, 
6001. It is obvious that the susceptibility barely changes 
with increasing lattice number in all the present data. 

From Figs. 5 and 6, one may also learn that the peak 
of reduced fidelity increases slightly as the system sizes L 
increases around the critical point J°=l/2. To study the 
scaling behavior of the fidelity susceptibility at the criti- 
cal point, we perform a finite-size scaling analysis. In Fig. 
7, dark lines with points indicate xf{Jz){x), the x— bond 
two-site reduced fidelity susceptibility at J z as a function 
of system size L, and the red circles show the results of 
fitting data. In the inset, the behavior of Xf(J%)(z), the 
z— bond two-site reduced fidelity susceptibility at J° is 
also shown. The reduced fidelity susceptibility at the 
critical point depends on the system size, which may be 
fit as 
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FIG. 7: (Color online) x-bond two-site reduced fidelity sus- 
ceptibility as a function of system sizes L at the critical point 
J^—l/2; dark lines with points indicate numerical data while 
open red circles indicate the fitting data. Inset: behavior of 
z-bond two-site reduced fidelity susceptibility as a function of 
the system size L at the critical point J^—l/2. 
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FIG. 8: (Color online) x-bond two-site reduced fidelity sus- 
ceptibility as a function of J x along the dashed line shown in 
the triangle for the various system sizes L=301, 1001, 3001, 
6001. 



As shown in Fig. 7, the fitting data agree with the orig- 
inal data rather well, so we may extrapolate Xf(Jz), 
in the thermodynamic limit. In the thermodynamic 
limit scales, and for the x— bond two-site reduced fidelity 
susceptibility, xf(J c z ) is 5.253(±0.001) , while xf{J z ) 
should be 16.430(±0.001) for the z— bond two-site re- 
duced fidelity susceptibility. 

These results confirm the previous conjecture; namely, 
that, the divergence of the global fidelity susceptibility is 
related to long-range correlations. The divergence disap- 



pears in the reduce fidelity susceptibility because the re- 
duced fidelity susceptibility averages over all long-range 
properties and only retains the nearest-site correlation 
information. 



IV. SUMMARY 



In summary, we study the reduced fidelity and reduced 

fidelity susceptibility in the Kitaev honeycomb model. 
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FIG. 9: (Color online) z-bond two-site reduced fidelity sus- 
ceptibility as a function of J x along the dashed line shown in 
the triangle for various system sizes L=301, 1001, 3001, 6001. 



It is shown that the nearest-two-site reduced fidelity 
susceptibility manifest itself as a peak at the quantum 
phase transition point, although the one-site reduced fi- 
delity susceptibility vanishes. Our results directly reveal 
that the reduced fidelity susceptibility is able to mark 
the quantum phase transition in the Kitaev honeycomb 
model, and thus suggest that the reduced fidelity suscep- 
tibility is still an accurate marker of the TPQT when it is 
properly chosen, despite its local nature. The conclusion 
that such a local quantity can characterize a TQPT is 
conceptually consistent with the fact that any physical 
observable is local in nature. 
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